Lifetimes of noisy repellors 



Holger Faisst and Bruno Eckhardt 
Fachbereich Physik, Phihpps Universitdt Marburg, D-35032 Marburg, Germany 

We study the effects of additive noise on the lifetimes of chaotic repellors. Using first order pertur- 
bation theory we argue that noise will increase the lifetime if the escape holes lie in regions where the 
unperturbed density is higher than in the immediate vicinity and that it decreases if the density is 
lower. Numerical experiments support the qualitative conclusions also beyond perturbation theory. 
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I. INTRODUCTION 



Noise can affect the behaviour of dynamical systems in 
many ways. It can induce transitions between otherwise 
disconnected regions (Kramers' theory, Q), it changes 
the scaling near bifurcations @|, it gives rise to stochas- 
tic resonance ^ and it can even change repellors into 
attractors 0. Some time ago Franaszek |^] studied the 
effects of additive noise on repellors |^ and found that in 
some cases it stabilized the dynamics, i.e. increased the 
life times. He studied this behaviour near crises and bi- 
furcations in the dynamics, but the reasons for the effects 
on the dynamics remained unclear. We here want to ap- 
proach the problem from the side of the attractor which 
then is perturbed to become a chaotic repellor. We will 
investigate the relation between the noise effects and a 
non-uniform density in the attractor. 

The hypothesis we want to test runs as follows: open- 
ing up the attractor into a repellor is achieved by punch- 
ing holes into the support of the attractor. Additive noise 
can push trajectories that would barely miss the holes 
into escape, but can also save trajectories that would es- 
cape in the unperturbed situation. Whether the life time 
increases or decreases then depends on which process is 
more likely: for a uniform density noise will kick out as 
many trajectories as it saves, so one cannot expect any 
effect. If the unperturbed density in the hole region is 
higher than in the immediate vicinity, more points will be 
saved than kicked out, and the lifetime should increase. 
If the unperturbed density in the hole region is lower, 
more trajectories will escape and the lifetime should be 
reduced. 

In section ^ we present the perturbative arguments, 
followed by numerical experiments in section lllll and 
some final remarks in section Hvl 



II. PERTURBATION THEORY 



We start with a d-dimensional map 



that has a chaotic attractor. The associated Frobenius- 
Perron equation for the evolution of densities p{x) is 



p„+i(x) ^ J dy(5(x-f(y))p„(y) 



p«(yi) 
l^f(y.)l 



(2) 
(3) 



where \D{\ is the Jacobi determinant and the summation 
extends over all points that map into x. Noise can be 
added as a Gaussian smearing of the propagator, as in a 
kicked system: the time evolution splits into two parts, 
the 'kick' i.e. the application of the deterministic 
map, and a free diffusive spreading. The free diffusion 
on the d-dimensional phase space is described by a diffu- 
sion kernel Ko that solves the appropriate free diffusion 
equation with (5-function initial conditions. For instance, 
for free diffusion in Euclidean space the Fokker-Planck 
equation for a density p is 



p = DAp 

and the diffusion kernel becomes 



-(y-x)V2Dt 



(4) 



(5) 



The combined evolution of kick and diffusion for some 
time T is then described by 

p„+i(x) J dzA'D(x,z,r) J dy5(z - f(y))p„(y|6) 

= J dyifD(x,f(y),r)p„(y) (7) 

dy/f(x,y)p„(y). (8) 



The evolution kernel if in (jHJ can be expanded in terms 
of left (A| and right eigenfunctions |A) with eigenvalues 
A, 



K = J2\\X){X\ 



(9) 



x„+i f(x„) 



(1) 



The existence of an invariant density on the attractor 
(the Sinai-Ruelle-Bowen measure) implies the presence 
of an eigenvalue A = 1 with right eigenvector the 



2 



invariant density, and a corresponding left eigenvector 
(1| = 1 because of conservation of probability. Since the 
kernel K is not self adjoint left- and right-eigenvectors are 
different. Simple examples show that the left eigenvectors 
develop fractal features 7]. With noise the finest scale 
structures are washed out, but higher levels of the fractal 
hierarchy survive. 

In order to turn the attractor into a chaotic repellor 
we punch holes into it. In applications the holes appear 
through crises and other perturbations and appear on 
many scales. For the purpose of the present analysis 
trajectories that enter the holes may be terminated, since 
their further evolution does not influence the escape rate. 
Let O be the domain over which trajectories are taken 
out and P the elimination projection: 



P(x,z) = 5(x- z) 



1 O 

p X e O 



(10) 
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This projection depends on a parameter p that will be 
useful in tests of the perturbation calculations: Trajec- 
tories entering O continue on with probability p and are 
taken out with probability \ — p. For the holes described 
before we have to take p = Q. The full evolution operator 
can then be written 

A>(x,y) =y"dzF(x,z)/v(z,y) (11) 
= i^(x,y)- / dz(J(x,z)-P(x,z))if(z,y)(12) 

JO 



(13) 



The last equation now has the form of an unperturbed 
propagator K plus a small perturbation aKi, where 
smallness is controlled by the localization in the region 
O and the rate 1—p with which points are taken out. a 
is a formal parameter that helps to organize the familiar 
perturbation expansion for eigenvalues and eigenvectors. 
The leading order result contains as usual the diagonal 
matrix element of the perturbation, 

A« Ao + a(A|A'i|A). (14) 

The deviations of the leading eigenvalue from Aq = 1 
define a decay rate, A = exp(— 7), where 

7«-a(A|ifi|A). (15) 

The systems we study here have piecewise smooth invari- 
ant densities and the first order correction to the decay 
rate has a regular dependence on the set O and the ex- 
traction rate p. In particular, for the logistic map with 
holes with p — that is studied in ^ expression p5|l 
gives the smooth background. The simulations by Paar 
and Pavin also indicate strong variations in escape 
rate near short periodic orbits which are connected to 
higher orders in perturbation theory and the complicated 
spatial structures that are characteristic for next to lead- 
ing eigenvectors (3] . 



FIG. 1: A graph of the map llfcil for c = 0.2 (dashed) and its 
invariant density p (continuous). 



III. PIECEWISE LINEAR MAPS WITH 
NON-CONSTANT INVARIANT DENSITIES 

Let us define a family of 1-d maps Xn+i ~ fd^n) with 
a parameter c S (0, 0.5), 



fcix) 



x/c < X < c 

l-2{x- c)c/(l - 2c) c<x <l/2 

c - (2a; - l)c/(l - 2c) 1/2 < x < 1 - c 

1 + (x - l)/c 1 - c < .X < 1 . 



(16) 

The invariant density is a solution to the Frobenius- 
Perron equation 



dyS[x - fc{v)]Pc{y) 



I f 



\^c{x^)\ ' 



(17) 
(18) 



where the sum is taken over all pre-images Xi of x. For 
(|16|l the normalized invariant density is 



Pcix) 



-r-rr — 7 < a; < c and 1 — c < a; < 1 

4c(l-c) (^g^ 



2(l-c) 



C < X < 1 — c . 



The main results do not depend on the specific value of 
c so that we can fix c = 0.2 and drop the subscript on / 
and p. 

For the eigenvalue and eigenvector analysis we use a 
matrix representation of the density evolution operator 
K . With the help of N characteristic functions 



(t>v{x) 



1, {v-l)/N<x<v/N,v^l,...,N 
0, elsewhere , 



(20) 
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FIG. 2: Eigenfunctions of map without noise and without 
holes. The upper row corresponds to the leading eigenvalue 
Ao = 1 and the lower one to the next to leading eigenvalue 
Ai = -0.6. 



FIG. 3: Escape rate 7 divided by the hole size e for the map 
llbll with holes but without noise. The dashed lines at 3.125 
and 1.125 mark the values expected within first order pertur- 
bation theory. 



B. With holes but without noise 



densities can be expanded as p = a,^(j),y (summation im- 
plied) and the evolution operator K becomes an x TV 
matrix. Densities are mapped according to 

= if.^a^r' , (21) 

and the matrix elements K^^ can be calculated from the 
images of the step functions: in the support of the char- 
acteristic function fi a uniformly distributed ensemble of 
5 X 10^ points is iterated and the probability to end up in 
the interval v is then the matrix element A'j/^ . The typi- 
cal size of matrices used in the calculations is iV = 3000. 

The lifetimes for the maps with holes can be obtained 
from the eigenvalues of the evolution operator after pro- 
jection onto the remaining intervals or directly from inte- 
grations of an ensemble of initial conditions. We followed 
10^ randomly selected initial points up to a maximal cut- 
off lifetime of 10'* iterations. The lifetime distributions 
decayed exponentially and the lifetimes estimated from 
this decay were within less than 1% of the eigenvalues. 



Pairs of holes of size e are added symmetrically at the 
edges of the invariant density. 'Outer' holes at (c — e, c) 
and (1 — c, 1 — c -|- e) lie within the high density region, 
'inner' holes at (c, c + e) and (1 — c — e, 1 — c) in the low 
density region. 

The perturbation theory from section^]predicts a lin- 
ear variation of escape rate with hole size, so that the 
ratio of escape rate to hole size should be constant. This 
is verified in Fig. |31 The numerical values for the escape 
rate 7 agree well with the expected values, sum of hole 
sizes times undisturbed density at holes, 7i^o = 2epi^o • 

With holes we no longer have conservation of proba- 
bility, and the left eigenvector to the leading eigenvalue 
will not be constant. It develops a fractal structure, that 
already for hole size e = 10"'^ is difficult to represent nu- 
merically. The left eigenstate for the leading eigenvalue 
Ao = 0.996855 for the map with outer holes is shown in 
Fig. 21 The largest eigenvalue corresponds to an escape 
rate 70 = 3.150 x 10"'^ in good agreement with the values 
extracted from Fig. 13 



A. No holes and no noise 

Without noise and without holes the map (|16|l maps 
the interval [0, 1] into itself. We expect one eigenvalue 1 
with right eigenvector the invariant density (as in Fig.[T]l. 
and a left eigenvector that is constant because of conser- 
vation of probability. The first two pairs of left and right 
eigenvectors are shown in Fig. |21 The next to leading 
eigenvalue is Ai = —0.6. Its left eigenvector shows the 
fractal structures one expects for such maps P]. 



C. Noisy map without holes 

For the map with noise we add Gaussian distributed 
random numbers at each time step, 



(22) 



where the ^„ are independent and identically distributed 
according to 



1 



(23) 
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FIG. 4: Left eigenfunction of map with outer holes of size 
e — 10~^ for the leading eigenvalue Aq. The hierarchy of 
peaks follows the pre-images of the holes. 



with (7 = \/2Dt the amplitude of the noise. Since it then 
is possible to leave the interval [0, 1] we close the system 
periodically by mapping points outside the interval back 
in using the modulo operation. The case with noise is 
the more regular one, and if the width of the character- 
istic function for the hole region is a fraction of the noise 
level the vectors converge rather reliably, as a compari- 
son between the results for matrix sizes TV — 3000 and 
N — 4000 for noise amplitude 10~^ showed. 

The eigenfunctions and eigenvalues in the presence of 
noise of amplitude a ~ 10^^ are shown in Fig. |31 The 
first eigenvalue remains at 1, the first left eigenvector 
is uniform, but the step in the first right eigenvector 
is smoothed out. With increasing noise amplitude this 
transition region becomes wider, as evidenced by the 
magnifications in Fig. 
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FIG. 5: Eigenfunctions for the noisy map without holes. The 
noise amplitude is a = 10~^. The upper row corresponds to 
the leading eigenvalue Ao = 1, the lower to the next to leading 
eigenvalue Ai — —0.6073. 




0.24 



D. Noisy map with holes 

We finally come to our model for a noisy repellor, the 
noisy map with holes. Fig. [7| shows for holes of size 
e — 1Q~^ the change in the escape rates as a function 
of noise amplitude. For outer holes the escape rate de- 
creases, for inner holes it increases, until for a noise am- 
plitude of about 10~^ they almost coincide. This is ex- 
actly what one would expect from the changes in invari- 
ant density shown in Fig. El the density on the upper level 
decreases and the one on the lower increases with the cor- 
responding changes in lifetime. If we had not closed the 
interval to a circle the loss of trajectories at the end of 
the intervals would have swamped this effect and the es- 
cape rate would have increased monotonically with noise 
level. 

The associated eigenfunctions are shown in Fig. |S1 for 



FIG. 6: Magnification of the invariant density for the noisy 
map without holes. Noise amplitudes are a — 10~^ (sharpest 
transition), 10""^, 10~^, and 10~^ (widest transition). 



outer holes of size 10~^ together with a noise amplitude 
lO^'^. Compared to the noise free case structures are 
smoothed out: for instance, the amplitudes of left eigen- 
functions decrease considerably. 



E. Failure of perturbation theory 

First order perturbation theory does not always work 
that well. Consider another map, Xn+i = g{xn) on the 
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FIG. 7: Escape rate 7 divided by the hole size e = IQ-^* for the 
noisy map with holes. When the holes are in the high density 
region (outer holes) the escape rate decreases, and when they 
are in the low density region (inner holes) it increases. When 
noise level a and hole size are comparable there is no diflerence 
between inner and outer hole placement anymore. 
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FIG. 8: Eigenfunctions of the map with noise amplitude a = 
10~^ and with outer holes of size e — 10~^. The eigenvalues 
are Ai = 0.99745 for the upper graph and A2 — —0.61138 for 
the lower one. 



interval [0, 1] with 

Its invariant density is (Fig. O: 
pix) 



_j 1/3 + 2a; < x < 1/3 
3(1 -a;)/2 1/3 < a; < 1 



3/4, < x < 1/3 
9/8, 1/3 < a; < 1 



(24) 



(25) 



Eigenfunctions for the two leading eigenvalues are 
shown in Fig. ^| Again we take a 'left' hole at [1/3 — 
e, 1/3] in the low density region and a 'right' hole at 
[1/3, 1/3 + e] in the high density region. First order per- 
turbation theory predicts the decay rates 7 as function 
of hole size e to be 



Ir 
11 



epr 
epi = 



: 9/8e 
3/4e. 



The numerical values extracted from Fig. 1111 are 



7r 



0.87e 
0.75e. 



(26) 
(27) 



(28) 
(29) 



In agreement with perturbation theory the escape rate is 
proportional to the hole size. The pre-factors from the 
two calculations agree for the left hole in the low, but 
disagree for the right one in the high density region. The 
deviation becomes smaller when the holes are partially 
closed, i.e., the parameter p in l(Tn|) is set to 0.99. Then 
the perturbative and numerical escape rates are in bet- 
ter agreement. This failure of the perturbative estimate 
seems to be closely connected with the presence of a pe- 
riodic orbit on the border of the interval. If the opening 
is shifted to slightly larger values then the slope agrees 
again with the perturbative results. 



IV. FINAL REMARKS 

Within the simple models studied here the hypothesis 
that the variation of lifetimes can be related to inho- 
mogeneities in the invariant density of the unperturbed 
attractor could be confirmed. Such inhomogeneities are 
most pronounced near bifurcations and crises 0, as in 
the work of Franaszek Q . The placement of holes from 
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FIG. 9; The map 1241 (dashed) and its invariant density (con- 
tinuous) . 
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the outside is less artifical than it may seem. In the 
case of riddUng bifurcations |0 Hne attractors are bro- 
ken up by holes that appear near periodic points that 
arc no longer transversally stable, and both position and 
widths can be controlled externally. Thus the observa- 
tions discussed here have some bearing on the effects of 
noise on the lifetimes in riddled attractors. 



FIG. 10: Eigenfunctions for the leading eigenvalue Ao = 1 and 
the next to leading eigenvalue Ai = —2/3 for map 12411 . 
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FIG. 11: Escape rates 7 divided by hole size e and by (1 — 
p) for the map 1241 1 vs. hole size for two placements of the 
hole and two extraction rates p. The dashed horizontal hnes 
mark the theoretical values 9/8 and 3/4. For the left hole the 
escape rate is in agreement with perturbation theory (circles, 
for p=0). For the opening to the right at p = (squares) the 
pre-factor does not agree with perturbation theory. For p = 
0.99 the values move up to the first order perturbation theory 
result (triangles). If the position of the opening is moved 
away from the critical point x = 1/3, e.g. to the interval 
[0.35, 0.35 + e], then the perturbative results is also obtained 
for p — (diamonds). 



